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ABSTRACT 

We report the discovery of KELT-7b, a transiting hot Jupiter with a mass of 1.28±0.18 Mj, radius of 
1.533^0.047 -ft-Ji and an orbital period of 2.7347749 ± 0.0000039 days. The bright host star (HD 33643; 

KELT-7) is an F-star with V = 8.54, T eS = 6789±fg K, [Fe/H]= 0.139^81, and lo Sff= 4.149 ±0.019. 

It has a mass of 1.535jig Q 54 M 0 , a radius of 1.732 ±q Q 4 and is the fifth most massive, fifth 
hottest, and the ninth brightest star known to host a transiting planet. It is also the brightest star 
around which KELT has discovered a transiting planet. Thus, KELT-7b is an ideal target for detailed 
characterization given its relatively low surface gravity, high equilibrium temperature, and bright host 
star. The rapid rotation of the star (73 ± 0.5 kms -1 ) results in a Rossiter-McLaughlin effect with an 
unusually large amplitude of several hundred ms -1 . We find that the orbit normal of the planet is 
likely to be well-aligned with the stellar spin axis, with a projected spin-orbit alignment of A = 9.7±5.2 
degrees. This is currently the second most rapidly rotating star to have a reflex signal (and thus mass 
determination) due to a planetary companion measured. 

Subject headings: planetary systems — stars: individual (KELT) techniques: spectroscopic, photo¬ 
metric 


1. INTRODUCTION 
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Transiting planets that orbit bright host stars are 
of great value to the exoplanet community. Bright 
host stars are ideal candidates for follow-up because the 
higher photon flux generally allows for a wider array 
of follow-up observations, more precise determination of 
physical parameters, and better ability to diagnose and 
control systematic errors. As a result, bright transit¬ 
ing systems have proven to be important laboratories for 
studying atmospheric properties of the planets through 
transmission and emission spectroscopy, for measuring 
the spin-orbit alignment of the planet orbi ts, a nd for de¬ 
termining precise stellar parameters f see I Win i (l20ll for 
a review). 

The Kilodeg ree Extremel y Little Telescope (KELT) 
transit survey (jPepner et all 1200 7) was designed to de¬ 
tect transiting planets around bright (8 < V < 10) stars. 
Very few (~ 3%) of the known transiting planet host 
stars are in this brightness range. This is because this 
range spans the gap between radial velocity surveys on 
the bright end, and the saturation limit of the majority 
of ground-based transit surveys on the faint end. The 
KELT-North (KELT-N) telescope targets this range us¬ 
ing a small-aperture (42 mm) camera with a a wide field 
of view of 26° x 26°. It observes 13 fields at declination of 
31.7 degrees, roughly equally spaced in right ascension, 
in total covering approximately 40% of the Northern sky. 
The KELT-N survey has been in operation since 2006 and 
candidates have been actively vetted since April 2011. 

The KELT-N survey has a lready announced four 
planet discoveries. KELT-lb (|Siverd et al.112012ft is a 
27 M.t brown dwarf trans iting a V = 10.7 F-star. KELT- 
2Ab dBeattv et all f20l2h is a hot Jupiter transiting the 
bright (V = 8.77) primary star in a visual binary sys- 
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Figure 1 . The top three panels show the effective temperature 
distr ibution for the star s targeted by the KELT-North transit sur¬ 
vey (jSiverd et al .1(201211 . the California Planet Survey (CPS) ra¬ 
dial velocity (RV) search (Wright et al. 2004), and the Kepler mis¬ 
sion (stars observed for all 16 quarters and with log g > 4.0 ac¬ 
cording to Kepler Q1-Q16 Stellar Parameters Database Q). The 
fourth and fifth panels from the top show the relative depth 
and mass of outer stellar convect ive zones at these temperatures 
(|van Saders Sz Pinsonneaultll200 ). while the si xth, bottom, panel 
show s the observed stellar usin i distribution ([Reiners & S chmitt! 
2003). The red dashed li ne at 6250K shows the approximate loca¬ 
tion of the Kraft Break dKraft!119700 . 

a http://exoplanetarchive.ipac.caltech.edu, NASA Exoplanet 

Archive 


tem. KELT-3b (jPepper et a,l.ll2013j) is a hot Jupiter tran¬ 
siti ng a V = 9.8 sligh tly evolved late F-star. KELT- 
6b dCollins et alJ[2014[ ) is a mildly-inflated Saturn-mass 
planet transiting a metal poor, slightly evolved late F- 
star. 

Because of its brighter magnitude range, the sample of 
host stars surveyed by KELT has a higher percentage of 
luminous stars than most transit surveys. This luminous 
subsample includes giants, as well as hot main-sequence 
stars and subgiants. Indeed, all five of the KELT-N dis¬ 
coveries to date (including KELT-7b) orbit F stars with 
T e g > 6100K. Such hot stars are typically avoided by 


radial velocity surveys. There is a transition between 
slowly and rapidly rotating stars known as the Kraft 
break ijKraftl fl970L 11967T) . Stars hotter than the Kraft 
break around T e g = 6250K typically have higher rota¬ 
tion velocities, making precision radial velocities more 
difficult. These higher rotation velocities reflect the an¬ 
gular momentum from formation, which is conserved as 
the stars evolve due to the lack of convective envelope. 
The lack of a convective envelope results in weak mag¬ 
netic fields and ine ff ective magnet i c brea king from stellar 
winds (e.g.. lvan Saders fe Pinsonneaultl ( 2013H L 

These points are illustrated in Figure [TJ which shows 
the distribution of effective temperatures for the KELT- 
N stellar sample, the sample of stars targeted by Kepler, 
and a representative radial velocity survey. The KELT- 
N targets plotted here are all of the bright (V < 11) 
putative dwarf stars in the survey which were selected 
by a reduced proper motion cut and with temperatures 
calculated from their J-K colors. Approximately 40, 000 
KELT-North targets (55%) are hotter than 6250K, 28 
of the CPS targets (2.3%) are hotter than 6250K, and 
approximately 20, 000 of the Kepler targets (20%) are 
hotter than 6250K. Also shown are theoretical estimates 
of the mass and radius of the convective envelope as a 
function of T e ,g for st a rs with solar metallicity and an 
age of 1 Gyr (|van Saders fe PinsonneauItll2012D . as well 
as the upper envelope of obs erved rotation velo c ities a s 
a function of T e ff based from IReiners fe Schmittl (1200311 . 

Hot stars pose both opportunities and challe ng es for 
transit surveys. On the one hand, hot starso have 
been relatively unexplored as compared to later spec¬ 
tral types. The first transiting pla net was discov- 
ered by radial velocity surveys dCharbonneau et aDl2000i ; 
iHenrv et all 12000 1. which initially targeted only late F, 
G, and iearly K stars. Due to the magnitude range of the 
stars surveyed and the choice of which candidates to fol¬ 
low up, the first dedicated ground-based transit surveys 
(lAlonso et al.ll2004lMcCullough et al.ll2006l : iBakos et al.1 

20071 iCollier Cameron et al.l l2007li ~ were also primarily 

sensitive to late F, G, and early K stars. As the value 
of transiting planets orbiting later stellar types was in¬ 
creasingly recognized, transit surveys began to survey 
lower-mass stars. Kepler extended the mag nitu de range 
of their target sample to fa inter magnitudes dGould et al.l 
l2003t IBatalha et alJlMflh in order to include a signifi¬ 
cant number of M dwarfs. The Kepler K2 mission will 
likely survey an even larger number of M stars than the 
prime mission (| Howell et al.l 120141) . MEarth is specifi¬ 
cally ta rgeting a sa mple of some 3,000 mid to late M 
dwarfs (jlrwin et al.1 120141) . Finally, HAT-South is sur¬ 
veying even fainter stars than HATNet, in order to in- 
c rease the fraction of late G, K, and even early M stars 
(|Bakos et al.H2013lf . 

As a result of this focus on later spectral types, the 
population of close-in, low-mass companions to hot stars 
is relatively poorly assayed. This is particularly true for 
stars which are both hot and massive; for example, only 
6 transiting planetary companions are known orbiting 


22 In this paper, we will follow IWinn et al.1 ( 20i7jj ) and define 
hot stars as those with T e g > 6250K. This is also roughly the 
temperature of the Kraft break for stars near the zero age main 
sequence (see [Kraft] 0.967;) and Figure [TJ. 
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stars with T e g > 6250K and M > 1 .5Mq 0 Build¬ 
ing a larger sample is particularly important given ex¬ 
isting claims that the population of planetary and sub- 
stellar companions to hot and/or massive stars is dif¬ 
ferent than that of cooler and less massive stars. In 
particular, there is evidence that hot Jupiters orbit¬ 
ing hot stars tend to have a large range of obliqui¬ 
ties (iWinn et al.ll201(H: lSchlaufmanll2010t lAlbrecht et al.1 
120121) . Based on surveys of giant stars, whose pro geni- 
tors are likely to be massive (I Johnson et al.1 (|2013l) , but 
see iLlovdl (|2013H b there have also been claims that the 
distribution of Jovian plan etary co mpani ons is a strong 
func tion of primary mass dBowler et al.1 (l2010f) . but see 
e.g. iMaldonado et al.1 (1201311 1. Finally, there is anecdo¬ 
tal evidence that massive substellar companio ns are m ore 
common around stars with T e g 6200K (iBouclrv et al.1 

[Mill . 

The challenges posed by hot stars are primarily due to 
the high rotation velocities. The large rotation veloci¬ 
ties of hot stars result in broad and weak lines, making 
precision radial velocity difficult. As a result, radial ve¬ 
locity surveys, and to a lesser extent transit surveys, have 
avoided targeting, or following up candidates from, such 
stars. Furthermore, for a fixed planet radius, the depths 
of planetary transits of hotter stars are shallower. This 
is exacerbated by the fact that stars with T e g > 6250K 
have lifetimes that are of order the age of the Galactic 
disk, and thus tend to be significantly evolved. 

However, there a number of ways in which these chal¬ 
lenges are mitigated for transit surveys. First, even 
though the transit depths are shallower, they are nev¬ 
ertheless greater than a few millimagnitudes, and thus 
readily detectable for Jovian-sized companions. There¬ 
fore, identifying such candidate transit signals is possi¬ 
ble even for main-sequence stars as hot as 7000K. Once a 
candidate signal is identified, its period can be confirmed 
with photometric follow-up. With a robust ephemeris in 
hand, radial velocity follow-up is greatly eased, as one is 
simply looking for a reflex variation with a specific period 
and phase (as opposed to searching over a wide range of 
these parameters, which increases the probability of false 
positives). Even with the relatively poor precision (a few 
100 ms -1 ) of radial velocity measurements of hot stars, 
it is possible to exclude stellar companions and detect the 
reflex motion of relatively massive planetary and substel¬ 
lar companions. 

Ultimately, however, it is precisely the high rotation 
velocities of hot stars that assist in robust confirmation 
of plan etary transits, via the Rossiter-M cLaughlin effect 
(RM) (iRossiterl 11924 iMcLaughlinl 11924) . The rotating 
host star allows one to measure the spectral aberration 
of the absorption lines due to the small blockage of light 
as the planet transits the rapidly rotating host. The 
magnitude of this effect can be directly predicted by the 
rotation velocity measured from the spectrum, combined 
with the transit depth and shape. The RM effect can 
therefore provide strong confirmation that the transit sig¬ 
nal is due to a planetary-sized object transiting the target 
star. However, for Jupiter-sized companions, this does 
not necessarily confirm the planetary nature of the oc- 
cultor, because low-mass stars , brow n d warfs , and Jovian 
planets all have roughly ~i?j (IChabrier & Baraffil2000( ). 


However, even a crude upper limit on the Doppler ampli¬ 
tude of a few kms -1 can then be used to exclude essen¬ 
tially all companions with masses in the stellar or brown 
dwarf regime. Thus, the Doppler upper limit, combined 
with the RM measurement, essentially confirms that the 
companion is a planet, i.e. that is both mass and radius 
are in the planetary regime. Furthermore, the shape of 
the RM signal allows one to measure the projected angle 
between the planet’s orbital axis and the star’s rotation 
axis. This projected ob liquity provides clues to the for¬ 
mation and evolution (|Albrecht et al.l 1 2012 ) history of 
hot Jupiters and substellar companions. This effect also 
provides an independent measurement of the rotational 
velocity of the star. 

In this paper, we describe the discovery and confirma¬ 
tion of a hot Jupiter transiting the bright V = 8.54 star 
HD 33643, which we designate as KELT-7b. In Section 
H we summarize the discovery photometric transit signal 
and the follow-up photometric and spectroscopic obser¬ 
vations. In Section [3l we discuss the analysis of the data 
obtained to determine stellar and planetary parameters. 
Section |4] considers the false positive scenarios and Sec¬ 
tion [5] discusses the results of this analysis. 


2. DISCOVERY AND FOLLOW-UP OBSERVATIONS 


2.1. KELT Observations and Photometry 

The KELT-N survey has a standard process of data 
reduction which will be bri efly describe d in this section. 
For more information see lSiverd et al.l (l2012h . KELT-7 
is in KELT-N survey field 04, which is centered on (a = 
05h: 54m: 14.71s, S = +31d: 39m: 55.10s; J2000). Field 
04 was monitored from 2006 October 26 to 2011 April 1 
collecting about 7800 images. The KELT-7 light curve in 
particular had 7745 points after a single round of iterative 
3 ct outlier clipping th at occurs just afte r the trend filter¬ 
ing algorthm (TFA) dKovacs et al.ll2005h . We reduced the 
raw survey data using a custom imple mentatio n of the 
ISIS image subtraction package (lAlard fc buntonl 119981 
IAlardll2f)OOD . combined wit h point-spread -function pho¬ 
tometry using DAOPHOT dStetson 119871). Using pr oper 
motions from the Tycho-2 catalog (iHog et af1l200C ) and 
J and H magnitu des from 2MASS dSkrutskie et al.1 20061 
iCutri et al.1 120031 ). we applied a reduced proper motion 
cut, (jGould et ah 120031) based on the i mplem entation of 
methods from iCollier Cameron et al.l (2007). This al¬ 
lowed us to select likely dwarf and subgiant stars within 
the field for further post-processing and analysis. We 
applied the TFA to each selected light curve to re¬ 
move systematic noise, followed by a search for tran¬ 
sit signals using the box-fit ting least squares algorithm 
(BLS) dKovacs et al.1 1200 2). For both TFA and BLS 
we used the versions found in the VARTOOLS package 
(lHartman et al.ir2008l) . 

One of the candidates from field 04 was star HD 33643 


/ 2MASS 05131092+3319054 / TYC 2393-852-1, located 
at (a = 05 h: 13to: 10.93s, 6 = +33d: 19m: 05.40s; J2000). 
The st ar has T ycho magnitudes Bg = 9.074 and Vt = 
8.612 (iHog et al.ll2000t ) and passed our initial selection 
cuts. The discovery light curve of KELT-7 is shown in 
Figure [2j We observed a transit-like feature at a period 
of 2.7347749 days, with a depth of about 8.28 mrnag. 


23 According to http://exoplanets.org 


2.2. Follow-up Time Series Photometry 
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Figure 2. Discovery light curve of KELT-7b from the KELT-N 
telescope. The light curve contains 7745 observations spanning 4.5 
years, phase folded to the orbital period of P = 2.7347749 days. 
The red line represents the same data binned at 1 hr intervals in 
phase. 



Time -T c (hrs) 


Figure 3. Follow-up transit photometry of KELT-7. The 
red over plotted line is the best fit transit model. The la¬ 
bels are as follows: MOR=University of Louisville Moore Ob¬ 
servatory; KEPCAM=KeplerCam at the Fred Lawrence Whip¬ 
ple Observatory; BOS=Bryne Observatory at Sedgwick (LCOGT); 
CROW=Canela’s Robotic Observatory; WHITIN=Whitin Obser¬ 
vatory at Wellesley College 


We obtained follow-up time-series photometry of 
KELT-7 to check for false positives and better determine 
the transit sh ape. We used the Tapir software package 
(j.Tensenll2013 f to predict transit events, and we obtained 
10 full or partial transits in multiple bands between Oc¬ 
tober 2012 and January 2014. All data were calibrated 
and processed using the AstroImageJ uackagtFI (AIJ; K. 
Collins ct al., in prep.) unless otherwise stated. 

We obtained one full transit of KELT-7b in the g-band 
on UT2012-10-04 at the University of Louisville’s Moore 
Observatory. We used the 0.6 m RC Optical Systems 
(RCOS) telescope with an Apogee U16M 4K x 4K CCD, 
giving a 26 x 26 field of view and 0.39 arcsec pixel -1 . 

We used KeplerCam on the 1.2 m telescope at the Fred 
Lawrence Whipple Observatory (FLWO) to observe two 
full 2 -band transits on UT2012-10-23 and on UT2012-11- 
03. We also observed a partial g-band transit on UT2012- 
11-22. On the night of UT2013-10-19 we observed a full 
i-band transit in combination with radial velocity ob¬ 
servations to measure the RM effect (described more in 
Section 12.31) . KeplerCam has a single 4K x 4K Fairchild 
CCD with 0.366 arcsec pixel -1 and a field of view of 
23.1 x 23.1 . The data were reduced using procedures 
outlined in lCarter et al.1 (120111) . which uses standard IDL 
routines. 

We observed a full transit in g-band on UT2012-11-14 
and a partial transit in Lband on UT2014-01-23 from 
the Byrne Observatory at Sedgwick (BOS), operated 
by Las Cumbres Observatory Global Telescope Network 
(LCOGT). BOS is a 0.8 nr RCOS telescope with a 3K 
x 2K SBIG STL-6303E detector. It has a 14.7' x 9.8' 
field of view and 0.572 arcsec pixel -1 . Data on the 
night of UT2012-11-14 were reduced and light curves 
were extracted us ing sta ndard IRAF/PyRAF routines 
as described in iFulton et all (|2011f ). Observations from 
UT2014-01-23 were analyzed using custom routines writ¬ 
ten in GOLF! 

We observed two full transits at Canela’s Robotic Ob¬ 
servatory (CROW) in Portugal. Observations were made 
using the 0.3 m LX200 telescope with a SBIG ST-8XME 
CCD. The field of view is 28 x 19 and 1.11 arcsec 
pixel -1 . Observations were taken on UT2012-12-08 and 
UT2013-01-29 in V-band and Lband, respectively. 

We observed one partial transit in /-band on the night 
of UT2013-01-27 at the Whitin Observatory at Wellesley 
College. The observatory uses a 0.6 nr Boiler and Chivens 
telescope with a DFM focal reducer that gives an effective 
focal ratio of f/9.6. The camera is an Apogee U230 2K 
x 2K with a 0.58 arcsec pixel -1 scale and a 20 x 20 
FOV. Reductions were carried out using standard IRAF 
packages, with photometry done in ALL 

Figure [3] shows each transit plotted with the best fit 
transit model over plotted in red. Figure [I] shows the 
combined and binned light curve with all 10 transits. 
Plots were generated during th e Global Fit analysis (See 
Section [Toll using EXOFAST (lEastman et al.ll2013f) . 

2.3. Spectroscopic Observations 

We used the Tilling hast Reflector Echelle Spectrograph 
f TRES: lFfires3l2008li . on the 1.5 m telescope at the Fred 


24 http://www.astro.louisville.edu/software/astroimagej 

25 GNU Data Language; http://gnudatalanguage.sourceforge.net/ 
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Figure 4. Top panel: All follow-up light curves from Figure 
[3] combined and binned in 5 minute intervals. This light curve 
is not used for analysis, but is shown in order to illustrate the 
best combined behavior of the light curve data set. The red curve 
shows the 10 transit models for each of the individual fits combined 
and binned in 5 minute intervals the same way as the data, with 
the model points connected. Bottom panel: The residuals of the 
binned light curve from the binned model in the top panel. 

Lawrence Whipple Observatory (FLWO) on Mt. Hop¬ 
kins, Arizona, to obtain spectra to test false positive sce¬ 
narios, characterize radial velocity (RV) variations and 
determine stellar parameters of the host star. We ob¬ 
tained a total of 64 TRES spectra between UT 2012 Jan¬ 
uary 31 and UT 2014 February 21. The exposure time 
varied from 90 — 2400 seconds depending on the weather 
conditions and the SNR we were trying to achieve. The 
spectra have a resolving power of R = 44 , 000 an d were 
extracted as described bv iBuchhave et all (l2010tl . 

Initially we obtained observations near phases 0.25 and 
0.75 in order to check for large velocity variations due to 
a small stellar companion responsible for the light curve. 
The spectrum appeared to be single-lined, and the veloc¬ 
ity variation that we saw was much too small to be due 
to a stellar companion so we continued observing to get a 
preliminary orbit. The orbit had a fair amount of scatter 
due to the rapid rotation of the host star so we stopped 
observing spectroscopically and opted instead to follow¬ 
up the star photometrically to confirm the depth, shape, 
and period of the transit, and to search for color-depth 
dependent depth variations indicative of a blended eclips¬ 
ing binary. Once we had observed multiple transits in 
different filters to determine that the transit depths were 
indeed achromatic, we started obtaining high signal-to- 
noise ratio spectra to refine the orbit and to determine 
stellar parameters of the host star. Table |T] lists all of 
the RV data for KELT-7. 

Of the 64 total spectra, 28 were taken on the night 
of UT 2013 October 19 to measure the RM effect and 
determine the projected obliquity of the system. Simul¬ 
taneous data were taken using the TRES spectrograph 
on the 1.5 m telescope and photometric data using Ke- 
plerCam on the 1.2 m telescope both atop Mt. Hopkins 
in AZ. We collected 28 RV spectra with a 9 minute expo¬ 
sure cadence and signal-to-noise ratio ranging from 127 
to 165 per resolution element. For the spectroscopic ob¬ 
servations we began collecting data an hour and a half 
before ingress but we only obtained 2 observations after 
egress due to morning twilight. Photometric observa¬ 
tions were gathered starting two hours prior to ingress 
until morning twilight which occurred about 10 minutes 
after egress. We obtained a total of 997 KeplerCam ob¬ 


servations at an exposure time of 2 seconds and a slight 
defocus of the image because of the brightness of the star. 

2.4. Adaptive Optics Observations 

We obtained adaptive optics (AO) imagery for KELT- 
7 on UT 2014 August 17 using the NIRC2 (instrument 
PI: Keith Matthews) with the Keck II Natural Guide 
Star (NGS) AO system (|Wizinowichl 1200(1) . We used 
the narrow camera setting with a plate scale of 10 mas 
pixel -1 . The setting provides a fine spatial sampling of 
the instrument point spread function (PSF). The observ¬ 
ing conditions were good, with seeing of 0.5". KELT-7 
was observed at an airmass of 1.31. We used a Br-y filter 
to acquire images with a 3-point dither method. At each 
dither position, we took an exposure of 0.5 second per 
coadd and 20 coadds. The total on-source integration 
time was 30 sec. 

The raw NIRC2 data were processed using standard 
techniques to replace bad pixels, flat-field, subtract ther¬ 
mal background, align and co-add frames. We did not 
find any nearby companions or background sources at 
the 5-a level (Fig. 0. We calculated the 5-cr detection 
limit as follows. We defined a series of concentric annuli 
centered on the star. For the concentric annuli, we cal¬ 
culated the median and the standard deviation of flux 
for pixels within these annuli. We used the value of five 
times the standard deviation above the median as the 5-cr 
detection limit. The 5-cr detection limits are Amag=2.5 
mag, 5.4 mag, 6.4 mag, and 7.3 mag for 0.1", 0.2", 0.5", 
and 1.0", respectively. 

3. ANALYSIS 
3.1. Stellar Parameters 

Using the Spectral Parameter Classification (SPC) 
(IBuchhave et al.ll2012ll technique, with T e g, log g, [m/H], 
and v sin i as free parameters, we obtained stellar pa¬ 
rameters of KELT-7 from the 64 TRES spectra. SPC 
cross correlates an observed spectrum against a grid of 
synt hetic spectra based on Kurucz atmospheric mod¬ 
els (llvuruczl 1 1992f ). The weighted average results are: 
T eff = 6779T50 K, log c?= 4.23T0.10, [m/H]= 0.12±0.08, 
and v sin i= 73.2 ±0.5 kms -1 . [m/H] was substituted for 
[Fe/H] in this analysis but we do not believe that this will 
affect the results. The weighted mean values were cal¬ 
culated by taking an average of the stellar parameters 
that were calculated for each spectra individually, and 
weights them accor ding to the cross-correlation function 
(CCF) peak height. 0 

3.2. Radial Velocity Analysis 

The relative RVs were derived by cross-correlating the 
spectra against the strongest observed spectrum from 
the wavelength range 4250 - 5650A. Figure [6] shows the 
best-fit orbit and computed bisectors with residuals. The 
best-fit orbit is a result of the EXOFAST Global Fit (See 
Section 13.51) assuming a fixed eccentricity of zero. The 
bisector analysis of the RVs taken out of transit showed 

26 SPC compares the observed spectra against a library of syn¬ 
thetic spectra calculated with the same mix of metals as the Sun. 
Since it uses all the lines in the observed spectra in the wavelength 
region covered by the library, the metallicity [m/H] is the same as 
[Fe/H] only if the mix of metals in the target star is the same as 
the Sun. 
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Table 1 

RV Observations of KELT-7 


Time 

BJDtdb 

Relative RV 
(ms -1 ) 

Relative RV 
(m s -1 ) 

error 

Bisector 

(ms -1 ) 

Bisector error 
(ms -1 ) 

Phase 

SNRe a 

2455957.836103 

102 


143 

-114 

209 

-0.3053 

68.8 

2455961.748540 

-472 


134 

-37 

143 

1.1253 

87.4 

2455963.607625 

-65 


144 

-35 

102 

1.8050 

122.7 

2455964.647614 

-387 


174 

170 

137 

2.1853 

82.2 

2455967.702756 

-129 


73 

16 

116 

3.3024 

101.0 

2455968.667152 

-25 


86 

-52 

61 

3.6551 

191.2 

2455970.667266 

-135 


136 

-79 

69 

4.3864 

106.3 

2455971.657526 

67 


95 

-199 

77 

4.7485 

161.3 

2455982.637552 

-13 


83 

127 

53 

8.7633 

153.3 

2455983.610040 

-150 


97 

55 

107 

9.1189 

174.2 

2455984.641133 

-197 


90 

86 

76 

9.4959 

180.2 

2455985.610681 

94 


104 

-55 

107 

9.8504 

122.0 

2455986.621278 

-267 


123 

-164 

156 

10.2200 

92.1 

2456019.663088 

-273 


115 

108 

92 

22.3017 

73.5 

2456020.642102 

-139 


107 

-100 

88 

22.6597 

125.5 

2456027.643080 

-52 


89 

-262 

119 

25.2196 

141.2 

2456310.790708 

55 


117 

144 

60 

128.7524 

296.9 

2456340.790829 

0 


48 

97 

44 

139.7219 

297.5 

2456584.819001 

-45 


88 

102 

131 

228.9507 

131.4 

2456584.826693 

81 


64 

26 

71 

228.9535 

127.1 

2456584.835357 

62 


100 

-42 

85 

228.9567 

146.6 

2456584.843054 

-1 


93 

75 

70 

228.9595 

148.1 

2456584.852250 

65 


85 

8 

95 

228.9629 

144.0 

2456584.859959 

-62 


63 

41 

87 

228.9657 

144.5 

2456584.867813 

13 


48 

14 

78 

228.9686 

145.7 

2456584.875811 

7 


68 

8 

104 

228.9715 

143.5 

2456584.884423 

96 


108 

-39 

74 

228.9746 

154.8 

2456584.892300 

230 


96 

-130 

57 

228.9775 

154.3 

2456584.899864 

181 


71 

-179 

70 

228.9803 

155.7 

2456584.907399 

265 


110 

33 

112 

228.9830 

161.2 

2456584.915097 

86 


82 

89 

61 

228.9858 

155.4 

2456584.923361 

146 


98 

138 

62 

228.9889 

154.1 

2456584.930798 

93 


81 

56 

43 

228.9916 

153.1 

2456584.938490 

-78 


65 

98 

91 

228.9944 

162.9 

2456584.946413 

1 


76 

-9 

55 

228.9973 

162.7 

2456584.956825 

-87 


103 

-80 

84 

229.0011 

161.8 

2456584.964968 

-172 


82 

-77 

78 

229.0041 

166.5 

2456584.972520 

-205 


109 

-70 

74 

229.0068 

165.2 

2456584.980362 

-255 


89 

-143 

82 

229.0097 

158.9 

2456584.987938 

-382 


80 

-23 

56 

229.0125 

164.0 

2456584.996897 

-610 


118 

132 

72 

229.0158 

162.3 

2456585.004641 

-437 


81 

240 

53 

229.0186 

165.4 

2456585.012118 

-223 


70 

336 

79 

229.0213 

162.9 

2456585.019665 

-86 


80 

207 

70 

229.0241 

159.9 

2456585.027143 

-37 


123 

243 

65 

229.0268 

156.1 

2456585.034609 

-226 


78 

145 

70 

229.0295 

156.2 

2456638.901881 

26 


98 

-169 

112 

248.7261 

198.1 

2456639.777025 

-123 


68 

-67 

73 

249.0461 

242.0 

2456640.913299 

-115 


80 

90 

53 

249.4616 

262.4 

2456641.702340 

-166 


78 

6 

67 

249.7501 

169.0 

2456642.742449 

-347 


66 

39 

92 

250.1304 

258.3 

2456693.611029 

72 


73 

-200 

72 

268.7305 

203.3 

2456694.641506 

-208 


89 

-67 

68 

269.1073 

186.2 

2456696.631120 

11 


73 

-149 

70 

269.8348 

229.6 

2456700.710497 

-519 


134 

-2 

75 

271.3264 

146.0 

2456701.756533 

-38 


67 

-206 

70 

271.7089 

238.9 

2456702.691941 

-207 


53 

-53 

50 

272.0509 

278.9 

2456703.670164 

-194 


59 

-35 

42 

272.4086 

263.7 

2456704.671919 

-94 


77 

61 

50 

272.7749 

219.1 

2456705.773764 

-316 


76 

-78 

75 

273.1778 

220.1 

2456706.652092 

-32 


68 

-41 

60 

273.4989 

254.3 

2456707.664631 

9 


63 

-63 

37 

273.8692 

238.4 

2456708.711942 

-381 


84 

80 

42 

274.2521 

206.1 

2456709.746023 

-21 


65 

-27 

59 

274.6302 

260.8 


a Signal to noise per resolution element (SNRe) which takes into account the resolution of the instrument. 
SNRe is calculated near the peak of the echelle order that includes the Mg b lines. 
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Table 2 

Stellar Properties of KELT-7 


Parameter 

Description 

Value 

Source 

Reference a 

Names 

O!J2000 


BD +33 977 

TYC 2393-852-1 
2MASS 05131092+3319054 
GSC 2393-00852 

HD 33643 

05 13 10.93 

Tycho-2 

i 

<5.72000 


+33 19 05.40 

Tycho-2 

i 

NUVqalex 


13.330 ± 0.905 

GALEX 

2 

Bjp 


9.074 ± 0.030 

Tycho-2 

1 

V T 


8.612 ± 0.030 

Tycho-2 

1 

V 


8.540 + 0.030 

SKY2000 

3 

B 


8.970 + 0.030 

SKY2000 

3 

U 


9.010 + 0.030 

SKY2000 

3 

Ic 


8.129 ± 0.051 

TASS 

4 

J 


7.739 ± 0.030 

2MASS 

5 

H 


7.580 ± 0.042 

2MASS 

5 

K 


7.543 ± 0.030 

2MASS 

5 

WISE1 


10.179 + 0.050 

WISE 

6 

WISE2 


10.844 ± 0.050 

WISE 

6 

WISE3 


12.766 + 0.180 

WISE 

6 

WISE4 


13.741 + 0.123 

WISE 

6 


Proper Motion in RA (mas yr" 1 ) 

10.40 ± 0.70 

UCAC4 

7 

V-5 

Proper Motion in DEC (mas yr" 1 ) 

-49.70 ± 0.60 

UCAC4 

7 

U h 

kms -1 

-33.5 + 0.2 

This paper 


V 

kms -1 

-9.7+ 1.8 

This paper 


w 

kms -1 

-8.4 ±0.9 

This paper 


d 

Distance (pc) 

129 + 8 

This paper 


A v 

Age (Gyr) 

Visual extinction 

1.3 ± 0.2 

0.13 ± 0.04 

This paper 
This paper 



a References: (l lH 0 g et al. (2000); (2 Martin et al. (2005); f3 lMvers e t al. (2001); (4)Richmond ct al. (2000); 
(5 lCutri et al.| (2003); Skrutsk ie et all ( 20061) ; (6 1 Wright et ahl (120101) ; ICutri & ct al. (2012); (7)Zacharias ct ah 
([ 2013 ) 

b 


Positive U is in the direction of the Galactic center. 






































Bieryla et al 


KELT-7 (HD 33643) N 

A 

E< - 


1 " 



Figure 5. Top: Br -7 AO image for KELT-7 (HD 33643) in Br -7 
(A = 2.1654 mm). North is up, east is left. The horizontal bar is 
1". No nearby companions or background sources were detected. 
Bottom: 5 a detection limit as a function of angular separation. 
Detecti on l imits at 0.1”, 0.2”, 0.5”, and 1.0” are also given in 
Section 12.41 

no indication that the bisector spans were in phase with 
the photometric ephemeris but the RMS was large due 
to the high v sin i. Despite the higher v sin i and resulting 
poorer precision, we were ultimately able to detect the 
reflex signal at high confidence (roughly 7a). We do see 
a correlation between the bisectors and RVs taken during 
the transit due to the RM effect. The relative RVs and 
bisector values are listed in Table [1] 

3.3. Rossiter-McLaughlin Analysis 

We performed an analysis of the RM data separately 
from the global fit analysis (see Section 1531) . To model 
the RM effect, we used parameter e stimation an d model 
fitting pro tocols as desc r ibed in iSanchis-Oieda et all 
(j2013l ) and 1 Albr echt et all (I2012|). The code implements 
formulas from IHirano et all (|2011l) , using the loss of light 
calculated from transit parameters and planet position as 
inputs. The transit data from the night of the RM event 
were used to determine the time of transit and transit 
parameters 6, R*/a, Rp/R+. Additional free parameters 
are v sin i and A to describe the amplitude and shape of 


400 

_ 200 

<0 

E 0 

> -200 

01 -400 
-600 

On 200 

°-§- -200 
200 

<n > 0 

m-S -200 

0.0 0.2 0.4 0.6 0.8 1.0 

Orbital Phase + 0.25 

Figure 6. TRES radial velocities of KELT-7. The green squares 
represent data taken on the night of the RM event, while the black 
circles are data that were not taken during transit. The phases 
have been shifted so that a phase of 0.25 corresponds to the time 
of the primary transit, Tq- Top panel: RV observations phased to 
the best orbital model with eccentricity fixed to zero and with no 
linear trend, shown in red. The predicted RM effect in the model 
shown incorporates the best fit model where A = 9.7 d= 5.2 degrees. 
Middle panel: Residuals of the RV observations to our circular 
orbital fit. Bottom panel: Bisector span of the RV observations as 
a function of phase. 

the signal, and a slope 7rm and offset 7rm to describe 
the orbital motion of the star. A is the angle on the 
sky measured clockwise from the sky-projection of the 
orbit angular momentum vector, to th e sky-projec tion of 
the stellar angular momentum vector (lOhta et al.ll2005h . 
The uncertainty of the model parameters were estimated 
using a Markov Chain Monte Carlo (MCMC) algorithm, 
where the number of chains was large enough to guaran¬ 
tee the robustness of the final values. 

The results from the analysis show that the sky- 
projected obliquity is A = 4. V 77 degrees. The analysis 
also gives us an independent measure of the projected 
rotational velocity, usm* = 66 (^g kins -1 , which is con¬ 
sistent with the SPC analysis (see Section 13.11) . The 
7rm offset was determined to be — 54+33 ms -1 . We 
can also use the result from the out of transit accel¬ 
eration 7 rm = —6711340 ms -1 day -1 to estimate the 
velocity semi-amplitude due to the planet. Using the or¬ 
bital period and assuming a circular orbit, we calculate 
Krv = 292 ± 146 ms -1 . The results from this analysis 
are shown in Figure [7] 

3.4. Time-series Spectral Line Profile 

The overall starlight that is blocked by the planet dur¬ 
ing transit will appear as a b ump in th e rota tional broad¬ 
ening function dCollier Cameron et al.ll2010D . Each spec¬ 
trum taken on the night of the RM event was cross- 
correlated against a non-rotating template. The CCFs 
of the out-of-transit data were median combined to cre¬ 
ate a master OOT CCF. The master was then subtracted 
from all of the in- and out-of-transit CCFs. Figure [ 8 ] is 
a greyscale plot showing the results. The bright white 
feature increasing in velocity with time is caused by the 
planet as it transits the star. This feature is what we 
would expect to see for a system that has low orbital 
obliquity: the planet crosses from the blue-shifted side 
of the stellar disk to the red-shifted side, and the center 
of the transit occurs at a r^sinmear zero. 

3.5. EXOFAST Global Fit 
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RM prediction for Kelt-7b 




0.85 0.90 0.95 1.00 1.05 1.10 

Time (BJD-2456584) 



Time - T c (hrs) 

Figure 9. RM results from the global EXOFAST fit with the re¬ 
flex velocity subtracted out. Top panel: RV data (green points) 
from UT 2013 October 19 with the best fit model shown in red. 
The two black points were taken on different nights and not in¬ 
cluded in the fit. The shape of the RM signal implies that the 
projected obliquity of the host star with respect to the planet is 
small. Bottom panel: The residuals of the data to the RM fit. 


Figure 7. RM results from our independent analysis. Spectro¬ 
scopic and photometric data are from UT2013 October 19. Top 
Panel: RV observations with the best RM fit shown in red. Bot¬ 
tom Panel: Photometric transit data with the best fit shown in 
red. Data are plotted in time for comparison. 



-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 

Vel. / Vsin(i) (70km/s) 

Figure 8. Time series of the residual average spectral line profile 
for data taken on UT2013 October 19. The bright white feature is 
caused by the planet transiting the star. 

We used a _ custom version of EXOFAST 

(jEastman et al.l 12013T) to determine a global fit of the 
system. EXOFAST does a simultaneous MCMC analysis 
of the photometric and spectroscopic data, including 
constraints on the stellar parame ters of M* a nd /A f rom 
the empirical Torres relations dTorres et al.l l2010h or 
Yons ei-Yale (YY) evolutionary models (iDemaraue et al.l 
12004 ). to derive system parameters. This method is 
similar to that described in detail in lSiverd et al.l ( 201 2b 
but we note a few differences below B 


27 In the EXOFAST analysis, which includes the modeling of 
filter-specific limb darkening parameters of the transit, we employ 
the transmission curves defined for the primed SDSS filters rather 
than the unprimed versions. We expect any differences due to that 
discrepency to be well below the precision of all our observations in 
this p aper and of the limb darkening tables from Claret & Bloemenl 

i2nm 


As initial inputs for EXOFAST we included as priors 
the orbital period P = 2.7347749 ± 0.000004 days from 
the KELT-N data and the host star effective temperature 
T e s = 6779 ± 50 K, metallicity [Fe/H] = 0.12 ± 0.08, and 
stellar surface gravity log g = 4.23 ± 0.10 from TRES 
spectroscopy. The priors were im plemented a s a x 2 
penalty in EXOFAST (see lEastman et ahl (l2013f) for de¬ 
tails). In fitting the TRES RVs independently to a Kep- 
lerian model we did not detect a significant slope in the 
RVs (i.e., due to an additional long-period companion), 
and we therefore did not include this as a free parameter 
in our final fits. 

We used the AIJ package to determine detrending pa¬ 
rameters for the light curves, such as corrections for air- 
mass and meridan flip. AIJ allows interactive detrending 
capabilities. Once we determined the detrending param¬ 
eters for each light curve, we fit the transit light curves 
using EXOFAST. All light curves are detrended by air- 
mass while the CROW light curve from UT2013 Jan¬ 
uary 29 was also detrended by meridan flip and average 
FWHM in the image. The raw data with detrending pa¬ 
rameters were used as the input for EXOFAST, and final 
detrending was done in EXOFAST. 

There were a few other considerations when running 
the global fit. First, we had to choose whether to include 
just the full transits, with an ingress and an egress, or 
to include all transits including the partial transits that 
were missing an ingress or an egress. Second, we had the 
option of allowing the orbital eccentricity and argument 
of periastron to float free or to fix them to zero and force 
a circular orbit. Third, we had to choose between con¬ 
straining the mass-radius relationship using the Torres 
relations or by using the Yonsei-Yale stellar models. Fi¬ 
nally, we had the option to include the RM observations 
and fit the RM RVs as part of the global fit. 

For the initial runs, we chose to use only the full tran¬ 
sits and the non-RM RVs to ensure convergence. We 
also chose to fix the eccentricity to zero and set the con¬ 
straint on the mass-radius relationship using the Torres 
relations. Once the fit converged, we ran it again but 
changed the constraint on the mass-radius relationship 
to the YY stellar models. We found the final parameters 
were in agreement within the uncertainties. This gave us 
confidence to start adding the partial transits. We ini- 
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Table 3 

Transit Times for KELT-7 


Epoch 

r P(~, 

( BJDtdb ) 

Error 

sec 

O-C 

sec 

O-C/Error 

Observatory 

-55 

2456204.817057 

64 

5.65 

0.09 

MOR 

-48 

2456223.959470 

31 

-83.91 

-2.70 

KEPCAM 

-44 

2456234.898861 

42 

-59.97 

-1.42 

KEPCAM 

-40 

2456245.839584 

50 

79.05 

1.56 

BOS 

-37 

2456254.045118 

63 

182.60 

2.88 

KEPCAM 

-31 

2456270.451621 

55 

-4.72 

-0.08 

CROW 

-13 

2456319.678871 

59 

102.16 

1.72 

WHITIN 

-12 

2456322.413721 

56 

108.34 

1.92 

CROW 

84 

2456584.950978 

47 

-19.45 

-0.41 

KEPCAM 

119 

2456680.667558 

87 

-77.12 

-0.88 

BOS 

F 


c 

"E" 

o 

I 

o 


BOS 

CROW 

KEPCAM 


O WHITIN 


6 
4 
2 
0 : 

-2 : 

-4 |. 

-150 -100 


-50 0 50 

Epoch 


100 150 


The residuals of the transit times from the best-fit 


observations at epoch -13 are hidden behind the CROW observa¬ 
tions at epoch -12. 


tially had trouble getting the fit to converge with all the 
light curves and found that we needed to add prior width 
constraint on the transit timing variation (TTV) and the 
baseline flux for the partial transits. We then released 
the prior on eccentricity and let it float free. The result 
was e = 0.013±g;gfg which was consistent with a circular 
orbit. 

The last step of the process was to include the 
RM velocities in the combined global fit. The RM 
data were modeled using the lOhta et ahl (1200511 an¬ 
alytic approximation. At each step in the Markov 
C hain, we interpolated the linear limb darkening tables 
of [Claret fc Bloemenl (2 01 ll : based on the chain’s value 
for log g, T e g , and [Fe/H] to derive the linear limb dark¬ 
ening coefficient, u. Our ground-based observations gen¬ 
erally do not constrain the two quadratic coefficients well 
enough to yield unique fits to both parameters. The un¬ 
certainties associated with this value are not true uncer¬ 
tainties because they include the covariances with the 
other parameters. For this reason we chose not to in¬ 
clude the linear limb darkening coefficient, it, in Table 0] 
The RM RVs were allowed a velocity offset (7 rm ) sep¬ 
arate from the non-RM RV dataset velocity zeropoint 
fay) . We allowed for this because stars h ave intrinsic 
jitter (j Albrecht et alJ 201 ^ : lWinn et al.ll2006H that can be 
significant in rapidly rotating stars F-stars. iCegla et all 
( 2014 ) have suggested that F-stars have been found to 
have more vigorous convective motions despite being 
magnetically inactive, and that the RV jitter is strongly 
correlated with the granulation flicker. We find that the 
offset between our 7 rm and "/rv values of ~ 100 m s ^ 1 is 
comparable to the RMS residual of the non-RM RV data. 
We also find the 7 rm offset determined in the global fit 
( 7 rm = — 35± —19 ms -1 ) to be consistent with the 7 rm 
offset determined in Section 13.31 The results from the 
EXOFAST RM fit are shown in Figure [9] 

For our final fits, we included all full and partial tran¬ 
sits, all the radial velocities including the RM velocities, 
and we assumed a circular orbit with no RV slope. The 
stellar and planetary values derived using the YY stellar 
models and using the Torres relation are shown in Table 
[4] for comparison. We chose to adopt the YY model as 
our fiducial values. We find that the spin-orbit alignment 
A = 9.7 ± 5.2 degrees, and our velocity semi-amplitude, 
Krv = 138 ± 19 ms -1 , determined from the global fit 
both agree with our independent solution discussed in 
Section [3~3l We also find that the resulting usini from 
the global fit (v sin i = 65+® kms -1 ) is in close agree¬ 


ment with the results from SPC and our independent 
RM analysis. 

The TTVs for all follow-up transits are shown in Fig¬ 
ure M The global fit Tc and P were constrained only 
by the RV data and the priors imposed from the KELT 
discovery data. Using the follow-up transit light curves 
to constrain the ephemeris in the global fit would artifi¬ 
cially reduce any observed TTV signal. As part of the 
global analysis, we fit as a free parameter a transit center 
time Tc for each transit shown in Table [3] A straight line 
was fit to all mid-transit times in Table [3j and shown in 
Figure 1101 to derive a separate ephemeris from only the 
transit data. We find T a = 2456355.229809 ± 0.000198, 
P = 2.7347785 ± 0.0000038, with a y 2 of 27.58 and 8 
degrees of freedom. While the y 2 is much larger than 
one might initally expect, this is likely due to systemat- 
ics in the transit data from the ground-based photome¬ 
try. Properly removing systematics in the partial transit 
data would be difficult, so we are therefore not convinced 
that this is evidence for TTVs. We were careful to check 
th at all t i mestamps w ere in BJDrdb time system using 
lEastman et al.l ( 2010 1 to convert timestamps. Further 
studies would be required to rule out TTVs. 

3.6. Evolutionary Analysis 

We use T e ff, logg, stellar mass, and metallicity derived 
from the EXOFAST global fits (see Section 13.51 and Ta¬ 
ble [4]) , in combination with the theoretical evo lutionary 
track s of the Yonsei-Yale stellar models (iDemaroue et al.l 
12004 . to estimate the age of the KELT-7 system. We 
have not directly applied a prior on the age, but rather 
have assumed uniform priors on [Fe/H], logg, and T e g, 
which translates into non-uniform priors on the age. Fig¬ 
ure EH shows the theoretical HR diagram (logg vs. T e g) 
with evolutionary tracks for masses corresponding to the 
± 1(7 extrema in estimated uncertainty. We adopt the 
Yonsei-Yale constrained global fit represented in the top 
panel. The estimated stellar mass (and secondarily the 
metallicity) define the model stellar evolutionary track 
from which the age is inferred in the HR diagram. Within 
the lcr uncertainties on the observed T e g , logg, and 
[Fe/H], the Yonsei-Yale evolutionary track gives an in¬ 
ferred age of 1.3 ± 0.2 Gyr. The bottom panel of Figure 
EH is shown as a c omparison using th e Torres constrained 
global fit values (|Torres et al.l 1201(1) . The Torres model 
provides empirical relationships between observed stellar 
parameters T e g, log< 7 , and [Fe/H], and stellar mass and 
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Table 4 

Median values and 68% confidence interval determined from EXOFAST global fit for KELT-7b 


Parameter 

Units 

YY-isocrone Values (adopted) 

Torres Values 

Stellar Parameters: 



M, 

Mass (Mq) 

1 roc + 0.066 

i.ooo _ 0 054 

1 yjqq + 0.069 

1.4oo_o 068 

R :» 

Radius (R©) 

i 709 + 0.043 

1 . 045 

1.715 ± 0.049 

L * 

Luminosity ( L© ) 

K 70 + O .37 
<o -0.36 

5 ni + 0,39 
°- o± -0.37 

P* 

Density (cgs) 

0 419+ 0 027 
u -^ xy — 0.025 

0 415+°- 029 

U.^tiO_o 026 

log g * 

Surface gravity (cgs) 

4.149 ±0.019 

4.140 ±0.019 

Trff 

Effective temperature (K) 

67891*9 

6789 ± 49 

[Fe/H] 

Metallicity 

0.13918;°™ 

nil q+0.080 
u.xxo_0 083 

v sin 7* a 

Rotational velocity (kms -1 ) 

65.0t°;g 

65-41*;® 

A 

Spin-orbit alignment (degrees) 

9.7± 5.2 

9.5l* : 2 

Planetary Parameters: 



p 

Period (days) 

2.7347749 ±0.0000039 

9 70477150 + 0 .0000040 

1 ' 1 ou -0.0000039 

a 

Semi-major axis (AU) 

0 0441 Cj+0-00062 

U. U4:4:10 _0 00052 

0 04S64“I" 0 00067 

U.U4:004:_o 00068 

M P 

Mass ( Mj) 

1.28 ±0.18 

1.25 ± 0.18 

R P 

Radius ( Rj ) 

1 roo+0.046 

j..ooo_o 047 

1 514+ 0 051 
1 - Oi4 _0.050 

PP 

Density (cgs) 

0 442 T 0 073 
u -^^-0.068 

0 446 +0 ' 074 

U.^o _ 0 069 

log gp 

Surface gravity 

9 101 +0.061 

o.-LO-L—o 068 

q 1 qn+0.060 

o.xou _0 068 

T eq 

Equilibrium temperature (K) 

2048 ± 27 

205112® 

e 

Safronov number 

0-04801“ 

0.0486 ± 0.0068 

©> 

Incident flux (10 9 erg s —1 cm -2 ) 

4,OOl8; 2 J 

4 02” 1 ” 0 ’ 22 
^.uz_o 21 

RV Parameters: 



T C 

Time of inferior conjunction (BJDtdb) 

2456223.9592 ± 0.0017 

2456223.9591 ±0.0017 

Krv 

RV semi-amplitude (ms -1 ) 

138 ± 19 

138 ± 19 

Krm 

RM semi-amplitude (ms -1 ) 

542l*J 

543l*i 

Mp sinz 

Minimum mass ( Mj ) 

1.28 ±0.18 

1.24 ±0.18 

Mp/Jlt. 

Mass ratio 

0.00080 ± 0.00011 

0.00081 ± 0.00011 

7 RV 

RV velocity zeropoint ms -1 

-133 ± 15 

-133 ± 15 

7 RM 

RM velocity zeropoint ms -1 

-35 ± 19 

-34 ± 19 

/(ml, m2) 

Mass function (Mj) 

0.00000080l°;°°°°°°i 

0. 000000801q; QQQQQQ29 

Primary Transit Parameters: 



Rp/R * 

Radius of the planet in stellar radii 

0 oq097T 0 00065 
u - uyub, ‘ -0.00064 

o.o9074i8;888?S 

a/R * 

Semi-major axis in stellar radii 

5 49+ 0 12 

o.^u_o.n 

5.47 ± 0.12 

i 

Inclination (degrees) 

83.761q3® 

83.72l®;« 

b 

Impact parameter 

0 597+° 0 22 

u.ou^ _ 0 025 

0 599 +0 023 

u.py »_ 0 026 

S 

Transit depth 

0.00828 ± 0.00012 

0.00823 ± 0.00012 

To 

Best-fit linear ephemeris from transits (BJDtdb) 

2456355.229809 ± 0.000198 

2456352.495016 ± 0.000191 

PTransit 

Best-fit linear ephemeris period from transits (days) 

2.7347785 ±0.0000038 

2.7347795 ± 0.0000037 

Tfwhm 

FWHM duration (days) 

0.12795 ±0.00046 

0.12821 ± 0.00047 

T 

Ingress/egress duration (days) 

n ni QQc+6.0u0y^ 

U.UIoijO—0.00089 

n m«zm+ 0 uuuy6 

U.Ul»4U_o 00093 

Tl4 

Total duration (days) 

0 14630” 1 ” 0 ' 00097 
u.moou _ 0 00092 

0.14662 ± 0.00098 

Pt 

A priori non-grazing transit probability 

0.1655 ±0.0034 

0 1662 +0 ' 0035 
u.iooz _ 0 0036 

Pt,G 

A priori transit probability 

0 1 QQ7+0-0043 

U.iyor-0.0042 

0 1993+ 0 0044 
u,iyyo _0.0045 

Secondary Eclipse Parameters: 



T S 

Time of eclipse (BJDtdb) 

2456222.5918 ± 0.0017 

2456222.5918 ± 0.0017 


a We adopted the SPC value for usinz (73 kms *) as our fiducial value since the EXOFAST RM analysis is not designed to model 
rapidly rotating stars. 


radius (see Table 12} . From this comparison we see that 
the age estimate quoted from the YY stellar model es¬ 
timate is consistent with that inferred from the Torres 
model estimated parameters to within ler. 

3.7. SED Analysis 

We construct an empirical spectral energy distribu¬ 
tion (SED) of KELT-7 shown in Figure IT21 We use 
t he n ear-UV bandpasses from GALEX ([Martin et al.1 
I2005D, the Bt and V t colors from the Tycho — 2 cat¬ 
alog (|Hpg et al.ll2000D . near-infrared (NIR) fluxes in the 


J a nd H passbands from the 2MASS Point S ource Cata¬ 
log dCutri et al.l 120031 ISkrutskie et al.l I2006D , and near- 
and mid-infrared fl uxes in the four WISE passbands 
(|Wright et al.M2010h to deriv e the SEP. We fit this S EP 
to the NextGen models from IHauschildt et al.l (1199911 by 
fixing the values of T e g, log g, and [Fe/H] inferred from 
the global fit t o the light curve and RV data as described 
in Section 13.51 and listed in Table 0] and then finding the 
values of the visual extinction Ay and distance d that 
minimize y 2 . We find Ay = 0.13 ± 0.04 and d = 129 ± 8 
pc with the best fit model having a reduced y 2 = 1.83. 





















12 


Bieryla et al 



T eff ro 



T eff ™ 


Figure 11. Theore tical HR diagrams bas ed on Yonsei-Yale stel¬ 
lar evolution models JDemarque et all20041) . The solid lines repre¬ 
sent the evolutionary tracks for the best-fit values of the mass and 
metallicity of the host star from the global fits using the Yonsei- 
Yale constraints (top panel) a nd th e Torres constraints (bottom 
panel) as described in Section 13.51 The tracks for the extreme 
range of la uncertainties on M* and [Fe/H] are shown as dashed 
lines, backeting the lcr range shown in grey. The red crosses show 
T e ff and log# from the EXOFAST global fit analysis, the blue 
dots represent the location of the star for various ages in Gyr. We 
adopt the Yonsei-Yale constrained global fit represented in the top 
panel resulting in an estimated age of 1.3 ±0.2 Gyr, where we note 
the uncertainty does not include possible systematic errors in the 
adopted evolutionary tracks. 



Figure 12. Measured and best-fit SED for KELT-7 from UV 
through mid-IR. The red error bars indicate measurements of the 
flux of KELT-7 in UV, optical, NIR, and mid-IR passbands and 
listed in Table [2l The vertical bars are the lcr photometric uncer¬ 
tainties, whereas the horizontal error bars are the effective widths 
of the passbands. The solid c urve is the best-fi t theor etical SED 
from the NextGen models of IHauschildt et al.l J19991) . assuming 
stellar parameters T e ff, log g, and [Fe/H] fixed at the adopted val¬ 
ues in Table l4l with Ay and d allowed to vary. The blue dots are 
the predicted passband-integrated fluxes of the best-fit theoretical 
SED corresponding to our observed photometric bands. T he 2 2 
micron band shows a slight IR excess as discussed in Section 13771 


The results from this analysis are shown in Table [2] We 
note that the quoted statistical uncertainties on Ay and 
d, are likely to be underestimated because we have not ac¬ 
counted for the uncertainties in values of T e gf, log < 7 , and 
[Fe/H] used to derive the model SED. Furthermore, it 
is likely that alternate model atmospheres would predict 
somewhat different SEDs and thus values of extinction 
and distance. 

Our SED analysis yields a slight IR excess i n th e 22 
micro n band which was also reported bvlMcDonal d et al.1 
( 2012 / in a study of IR excess of Hipparc os st ars. Due 
to the young age of this star (see Section Id.61) , the de¬ 
tection of this excess could be evidence for a debris disk, 
though we suspect that it is likely due to background 
nebulosity. Inspection of the WISE image in the 22 mi¬ 
cron band shows clear background nebulosity associated 
with a nearby bright embedded star forming region that 
is very bright in teh WISE 22 micron image. Therefore 
we consider it likely that this background nebulosity is 
the cause of the apparently slight excess in the WISE 22 
micron passband. In any event, the excess is only ~ 2a 
and appears only in this one band, therefore it has no 
impact on the overall SED model fit. 

3.8. UVW Space Motion 

We evaluate the motion of KELT-7 through the Galaxy 
to place it among standard stellar populations. The ab¬ 
solute heliocentric radial velocity is +39.4 ± 0.1 kms -1 , 
where the uncertainty is due to the systematic uncer¬ 
tainties in the absolute velocities of the RV standard 
stars. Combining the absolute TRES RV with the 
distance from the spectral energy distribution analysis 
and proper motion inform ation from the UCAC4 cata¬ 
log (|Zacharias et al.l|2013f) . we find that KELT-7 has a 
U, V, W (where positive U is the direction of the Galactic 
center) of —33.5 ± 0.2, —9.7 ± 1.8, —8.4± 0.9, all in units 
of kms -1 , making this a thin disk star dBensbv et al.l 
l2003lb 

4. FALSE POSITIVE ANALYSIS 

There are many signals that could be mistaken for a 
planetary transit, so it is important to address some of 
these false positive scenarios. There are several reasons 
to favor a planetary signal over a false positive scenario 
for KELT-7b. 

KELT has a very small aperture, and thus a very large 
point spread function (PSF), so many initial detections 
turn out to be blended starlight from more than one star 
in the PSF mimicking a transit signal. Therefore, it is 
important that we follow up our initial detection with 
seeing-limited telescopes (i.e., with PSFs of ~ 1") to 
rule out any blended eclipsing binaries. Observations us¬ 
ing larger telescope in multiple filters then resolve stars 
that are blended even at the 1 " resolution, which typi¬ 
cally turn out to be bound systems such as hierarchical 
triples. Our follow-up transits were observed in several 
different bandpasses ( Vgriz ), and we found no evidence 
of a wavelength-dependent transit depth. 

We carefully inspected our spectra to look for light 
from another source. We did not see any evidence for 
the spectrum being double- or triple-lined. Our bisector 
analysis of the RVs taken out of transit showed no indi¬ 
cation of being in phase with the orbital solution but we 
do see a correlation between bisector variation and RV 
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Figure 13. Transit depth as a function of apparent V magnitude 
of the host star for a sample of transiting systems with bright 
(V < 11) hosts. KELT-7b is shown as a pink diamond. Bright 
stars with deep transits are generally the best targets for detailed 
follow-up. 

variation of the spectra taken during transit due to the 
RM effect. 

Our global fit with all spectroscopic and photometric 
data is well modeled by that of a transiting planet around 
a single star. We find that the log g derived from our 
global fit, 4.149±0.019, is consistent within errors to log g 
derived from our SPC analysis, 4.23±0.1. The amplitude 
of the RM signal is consistent with that expected from 
the v sin i measured from the stellar spectrum and the 
depth and impact parameter measured from the high- 
precision transit light curves (see Section 1531 and Section 

1531) . 

Finally, we obtained AO images, which exclude com¬ 
panion sources beyond a distance of 0.1", 0.2", 0.5" and 
1.0" from KELT-7 down to a magnitude difference of 2.5 
mag, 5.4 mag, 6.4 mag and 7.3 mag respectively, at a 
confidence level of 5a (see Figure [5]). 

We conclude that all the evidence is best described by 
a transiting hot Jupiter planet orbiting a rapidly rotating 
F-star. There is no significant evidence suggesting that 
the signal is better described from blended sources. 

5. DISCUSSION 

We have presented the discovery of KELT-7b, a hot 
Jupiter planet orbiting the ninth brightest star to host a 
known transiting planet. This is the fifth most massive 
star and fifth hottest star to host a transiting planetPT 
Figure[l3]shows the V magnitude versus transit depth for 
known transiting systems with V< 11, with the KELT 
discoveries highlighted. KELT-7b is an excellent candi¬ 
date for future detailed atmospheric studies because it 
is a bright host star and it has a relatively deep transit. 
Although we suspect that it may be due to background 
nebulosity because the star lies in a region of consid¬ 
erably higher IR nebulosity than its surroundings, the 
slight IR excess we find at 22 microns can be confirmed 
or excluded using follow up observations. 

28 According to http://exoplanets.org/ 


KELT-7 is a hot (~ 6800 K), rapidly-rotating (vsini 
~ 73 kms^ 1 ) star, and its planetary companion was 
originally confirmed via the RM effect, which was eas¬ 
ily detected with an amplitude of several hundred ms -1 . 
On the other hand, the reflex radial velocity motion of 
the star due to the companion was much more difficult 
to detect, although we did ultimately detect the signal 
at high confidence and to date this is the second most 
rapidly rotating transiting system to have this motion 
measured. This discovery therefore illustrates both the 
opportunities and challenges associated with confirming 
planetary companions transiting hot stars. We note that, 
because of its brighter magnitude ra nge, the Tran siting 
Exoplanet Survey Satellite (TESS) (iRicker et all 1201411 
will also survey a large number of hot stars. Therefore, 
at least some of the lessons learned from KELT for char¬ 
acterizing the population of planets orbiting hot stars are 
likely to apply to TESS as well. 

The fast rotation also allowed us to measure the RM 
effect and measure the projected obliquity of the sys¬ 
tem. Understanding the projected spin-orbit alignment 
of a planet and host star can allow us to infer informa- 
tion about the for matio n and evoluti o n of h ot Jupiters. 

I Winn et~aTI (120101) and iSchlaufmanl (|2010l) . by differ¬ 
ent methods, proposed that hot stars ( T e g > 6250 K) 
that host a tran siting hot J upiter typically have high 
stellar obliquity. iWinn et~ahl (|2010h suggested that hot 
Jupiter systems initially have a broad range of obliqui¬ 
ties, but the cool stars eventually realign with the orbits 
of their companions because the y undergo mor e rap id 
tidal dissipation than hot stars. 1 Albrecht et al.l (120121 ) 
did an R M analysis on a sample that nearly doubled 
the lWinn et al.l (120101) sample and confirmed the correla¬ 
tion of projected o bliquity a nd the effective temperature 
of the star. lAlbrecht et al.l (120121) also showed that the 
obliquity of systems with close-in massive planets have a 
dependence on the mass ratio and the distance between 
the star and planet. Specifically, they found that higher 
obliquities are measured in systems where the planet is 
relatively small. 

The KELT-7 system consists of a transiting hot Jupiter 
on a fairly close orbit (a = 0.04 AU) to its massive and 
hot host star. We measured the system to have a low 
stellar obliquity (A = 9.7 ± 5.2 degrees). One might ex¬ 
pect that this planet formed with a low obliquity and mi¬ 
grated in close to the star because it has been suggested 
that if the planet formed around a hot host star with a 
high obliquity it would be unable to realign due to the 
lack of convective envelope. With a larger sample of hot 
stars with transiting planets with projected obliquities, 
it will become possible to disentangle the dependences 
of stellar effective temperature, age, planet mass, and 
orbital distance on the projected obliquity. Ultimately, 
this will enable a deeper understanding of how systems 
form and evolve over time, and allow us to distinguish 
which systems are truly unique. 
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